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Abstract 

We prove positivity results about linearization and connection coeffi- 
cients for Bessel polynomials. The proof is based on a recursion formula 
and explicit formulas for the coefficients in special cases. The result implies 
that the distribution of a convex combination of independent Student- 
t random variables with arbitrary odd degrees of freedom has a density 
which is a convex combination of certain Student-t densities with odd 
degrees of freedom. 
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1 Introduction 

In this paper we consider the Bessel polynomials q n of degree n 

n 

q n (u) = Y J ^ ) u k , (1) 

where 



(2) 



(n) = (t) 2* = n\{2n-k)\2 k 
° k {2n)\(n-k)\kV 
The first examples of these polynomials are 

u 2 

q Q (u) = 1, qi (u) = 1 + u, q 2 (u) = 1 + u + y . 
They are normalized according to 

Qn (0) = 1, 

and thus differ from the polynomials 9 n (u) in 9 by the constant factor 
i.e. 

(2nV 

°n(u) = ——q n (u). 
n\l n 

For v > we recall that the probability density on R 



f v {x) = v -^, A„= "Z?[ (3) 
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has the characteristic function 

e ix yf u (x)dx = K(\y\), y G M, (4) 

J — oo 

where 

k u (u) = Yj^jU u K u (u) , u>0, (5) 

and K v is the modified Bessel function of the third kind. If v = n + \ with 
n = 0, 1, 2, . . . then 

fc„(u) = e-"g n (u), u > 0, (6) 

and is called a Student-t density with 2v = 2n + 1 degrees of freedom. For 
v = \ then is density of a Cauchy distribution. Note that for simplicity we 
have avoided the usual scaling of the Student-t distribution. 

In this paper, we provide the solutions of the three following problems: 
assuming in the rest of the paper that a is a constant with < a < 1, 

1. positivity and explicit values of the connection coefficients cl (o) in the 



expansion 

n 

(n) 



In (au) = E c k (°) 1^ {u) , (7) 



2. positivity and explicit value of the linearization coefficients j3\ (a) in the 
expansion 

n 
i=0 



q n (au)q n ((l - a)u) = E I3\ n \a)q n+i (u) , (8) 



3. positivity of the linearization coefficients P^ l ' m \a) in the expansion 

n+m 

q n (au)q m ((l - a)u) = E /3^' m \a)q k (u) . (9) 

k=nAm 

Note that /3 4 (n) (a) = /?i+?( a ) and that © 

is a special case of @ corresponding 
to m = with 4 n) (a) = ^ n,0) (o). Note also that n = in © yields 



E /f ,m) («) = i, 

fe=nAm 



so is a convex combination. As polynomial identities, ©-© of course hold 
for all complex a, u, but as we will see later, the positivity of the coefficients 
holds only for < a < 1. 

Because of @ and (JHJ) formula © is equivalent with the following identity 
between Student-t densities 

1 i / \ n+m 

~f 4.1 P) *T L ~f , 1 (t— ) = E /t' m) (10) 

a «+2 va/ I - a ™+2 Vl - a/ ^— ' k +2 

k=nAm 
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for < a < 1 and * is the ordinary convolution of densities. 

Although @ is more general than (JTJ),©, we stress that we give explicit 
formulas below for c^\a) and fl^ (a) from which the positivity is clear. The 

positivity of /3^. n ' m ^ (a) for the general case can be deduced from the special cases 
via a recursion formula, see Lemma 13.41 below. 

These problems have an important application in statistics: the Behrens- 
Fisher problem consists in testing the equality of the means of two normal 
populations. Fisher [I] 1 has shown that this test can be performed using the 
d— statistics defined as 

df 1 j 2 £ = t\ sin6> — t 2 cos 8, 

where t\ and i 2 are two independent Student-t random variables with respec- 
tive degrees of freedom f\ and / 2 and 9 £ [0, |] . Many different results have 
been obtained on the behaviour of the d— statistics. Tables of the distribution 
of df lt f 2t g have been provided in 1938 by Sukhatme at Fisher's suggestion. 
In 1956, Fisher and Healy explicited the distribution of df x f 2 $ as a mixture of 
Student-t distributions (Student-t distribution with a random, discrete num- 
ber of degrees of freedom) for small, odd values of /i and / 2 . This work was 
extended by Walker and Saw who provided, still in the case of odd num- 
bers of degrees of freedom, an explicit way of computing the coefficients of the 
Student-t mixture as solutions of a linear system; however, they did not prove 
the positivity of these coefficients, claiming only 

"Extensive numerical investigation indicates also that rji > for 
all i; however, an analytic proof has not been found." 

This conjecture is proved in Theorem 2 and 3 below. Section 2 of this paper 
gives the explicit solutions to problems 1, 2 and 3, whereas section 3 is dedicated 
to their proofs. The last section gives an extension of Theorem 2 in terms of 
inverse Gamma distributions. 

Using the fact that the Student-t distribution is a scale mixture of normal 
distributions by an inverse gamma distribution our positivity result is equivalent 
to an analogous positivity result for inverse gamma distributions. This result 
has been observed for small values of the degrees of freedom in ^3]. In jS] the 
coefficients are claimed to be non-negative but the paper does not contain any 
arguments to prove it. 



2 Results 

2.1 Solution of problem 1 and a stochastic interpretation 
Theorem 2.1 The coefficients (a) in write 

(n-fc)A(fc+l) , , N 7 , ,\ 

4"'c>^§ e G;tl r )( B ;*r>-)' 



1 the collected papers of R.A. Fisher are available at the following address 
|http : / / www-library, adelaide.edu.au / digitised /fis her / 1 
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for < k < n — 1 while c^ \a) = a n and, as < a < 1, i/iey are positive. 

A stochastic interpretation of Theorem 12 . II writes as follows: replacing u by \u\ 
and multiplying equation by exp(— |u|), we obtain 

n 

e -^a)\u\ e -a\u\ qn (a|u|) = C W (a) g fc (H)e-l"l. (11) 

fc=0 

Equality (jllj) can be interpreted as follows: the convex combination of an in- 
dependent Cauchy variable C and a Student-t variable X n with 2n + 1 degrees 
of freedom follows a Student-t distribution with random number 2K (to) + 1 of 
degrees of freedom: 

(1 - a) C + aX n = X K ^) , 
where K (w) E [0, n] is a discrete random variable such that 

Pv{K (w) = k} = cj. n) (a), < k <n. 

2.2 Solution of problem 2 and a probabilistic interpretation 
Theorem 2.2 T/ie coefficients /3- (a) in |0) write 

= (4a(l- ))-f 7 ^ T V2- 2 " <2n - 2i)!(2 " + 2i)! 



(2n)!y (re-i)!(n + i)! 



n—i 



and, as < a < 1, i/iey are positive. 

A probabilistic interpretation of this result writes as follows. 

Corollary 2.3 With a = s'mO, f\ = f% = 2n + 1, statistic df lt f 2t g follows 
a Student-t distribution with a random number of degrees of freedom F (lo) 
distributed according to 

Pr {F (u) = 2n + 2i + 1} = (3\ n) (a), < i < n. 
2.3 Problem 3 

Theorem 2.4 The coefficients ^' m \a) in ^ are positive for < a < 1. 

We were unable to derive the explicit values of the coefficients (a): how- 

ever, their positivity allows us to claim the following Corollary. 

Corollary 2.5 With a = sin#, f\ = 2n + l,/2 = 2m + 1, statistic df x j 2 ^ 
follows a Student-t distribution with a random number of degrees of freedom 
F (uj) distributed according to 

Pr {F (cj) = 2k + 1} = /#' m) (a) , nAm<k<n + m. 
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Theorem 2.6 For k > 2 let n\, . . . ,n/~ be nonnegative integers and let a%, . . . , 

be positive real numbers with sum 1. Then 

L 

Qui {0-lU)qn2 

(a 2 u) ■ ■ ■ q nk (a k u) = '^2f3 j q j (u), u€R (12) 
i=l 

with nonnegative coefficients (3j with sum 1 and I = min(ni, . . . , n k ), L = n\ + 
h n k . 



3 Proofs 

3.1 Generalities about Bessel polynomials 

As a preparation to the proofs we give some recursion formulas for q n . They 
follow from corresponding formulas for 9 n from [§], but they can also be proved 
directly from the definitions (^Q) and ((21) • The formulas are 

u 2 

gn+l(tt) = 1n(u) + 2 _ Qn-^U), U > 1, (13) 

XL 

qM = q n (u) ~ 2n _ 1 gn-i(it), n>l. (14) 



We can write 

n 

u n = Y^8<( l) q t {u), n = 0,l,... (15) 

i=0 

and 5} is given by a formula due to Carlitz see [01 p. 73] or [TT] : 

f (n+l)! (-1)"-'(2Q! for n^l < - <n 
_ J 2 " (n-i)W!(2i+l-n)! 101 2 - 1 - n (16) 
\ for < i < 2fi 

Later we need the following extension of (|13|) which we formulate using the 
Pochhammer symbol (z) n := z(z + 1) • • • (z + n — 1) for 2; £ C, n = 0, 1, . . . . 

Lemma 3.1 For < k < n we have 

k 

2k 1 \ ( n ' fc ) I \ 

i=0 

where 

Jn,k) _ r,2k (k\ fm h , 



7 p) = 2 2 fe (fKn-fc + iJ^C-n -i)*-,. (17) 



Proof: The Lemma is trivial for = and reduces to the recursion (|13|) for 
k = 1 written as 

n 2 g n _i(ii) = 2 2 (n - |) 2 (<fo+l(«) - • (18) 
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We will prove the formula (|17|) by induction in n, so assume it holds for some 
n and all < k < n. Multiplying the formula of the lemma by u 2 we get 



u 



2k+2 



q n ~k{u) =^2% n ' k \ 2 q n +i(.u), 



i=0 



hence by lfT5)l 



,2(*H-1), 



9n+i-(fc+i)M = 22 ( n + i + |)a fen+i+2(n) - g n+i+ i(u)] 



i=0 



= ^' k) 2 2 (n + k + \) 2 q n+k +2{u) 
k 

+ £ 2 2 (n + i + |) [ 7 £f >(« + i - |) - 7 f ' fc V + <+§)] 
i=l 

- 7 ^2 2 (n+|) 2 g n+1 ( U ). 
Using the induction hypothesis we easily get 

(n,fc) 2/ , ; , 1\ n 2k+2/ 7 1\ (n+l,fc+l) 

and 

- 7 ^2 2 (n + §)(n + i) = 2 2fc + 2 (n - k + ±) fe+1 (-n - |) fc+1 = . 
Concerning the coefficient C to q n+ \ + i(u) above we have 

C = 2 2k+2 (n + i + { k \{n-k + \) k+i _ x (-n - h)k-i+l(n + 



(n-k + \) k+i (-n - \)k~i{n + «'+§) 

V * / 

2k+2 {n-k + \) k+1+i {-n-\) k . 



(— n — \ + k — i) 



2/k-i 

k 



i - 1 

2 2 *+ 2 ( n -k+ i) fc+1+i (_„_ i) fc 

2 2 fc+ 2^+l^ (n _ A;+ i w( _ 



3\ (n+l,fc+l) 



□ 



We stress that Lemma IM.ll is the special case v = n + | of the following 
recursion for modified Bessel functions of the third kind. 

Lemma 3.2 For all v > and a// nonnegative integers j < v we have for 
u > 



'^K v ^ (u) = £>2)^ 



i=0 



r(!/ + i) 



r(i/ + 1 - 0" - *)) 



u^i^ +i (u) 
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and 



^ Vv r (^ + 1 - u - «)) r (^- j) 

Proof: The second formula follows from the first using formula (J5J), and 
the first can be proved by induction using the following recursion formula for 
modified Bessel functions of the third kind, cf. \Y2i p. 79] 

K v -x(u) = K v+1 (u) - —K v (u). 

u 

We skip the details. □ 

3.2 Proof of Theorem EH] 
From (UJ and (JT3J) we get 

n j n 

Qn {au) = ^2 a f )a3 Yl 5 i qi ^ = Yl 4 ( a )?* ( u ) 
j=0 i=0 k=0 



with 



n 
j=k 



fc ^!_(2fc)! ^ _ fc ( 2ra -j)!(j + l) 

(2")l k\ ,_J-?^> ] (n-j)\(j-k)\(2k + l-j)\ 



In particular c$ (a) = a n and for < k < n — 1 



where 



(n-fe)A(fc+l) . . 

/ si (2n — « — «)!(« + z + 1) 



(n-k- i)\i\(k + 1 - i)\ 



We clearly have 



(„-fc)A( fc+1 ) (r) 

P (a)= (-ir^pa-af 

with 

(n-fc)A(fc+l) . 

[n — k — — ry\k A- I — i)\ 

and we only consider < r < (n — fc) A (A; + 1). To sum this we shift the 
summation by r. For simplicity we define T := (n — k — r) A (A; + 1 — r) and get 

(-lYv^(l) = Vf IV {2n-k-r-i)\{k + r + i + l) 
L) p W L > ( n _ k-r-i)\i\{k + l-r-i)V 
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We write k + r + 1 + i = (2k + 2) — (k + 1 — r — i) and split the above sum 
accordingly 

(-l)V r) (U = (2* + 2) £(-!)' (2n-ib-r-<)! 



(n — A: — r — i)\ i\ {k + 1 — r — i)\ 



i=0 

(2n — — r — i)\ 



(n — k — r — (k — r — i)\ 
Note that for nonnegative integers a, b, c with b, c < a we have 

Z^" 1 ) ( 6 _i)!( c _i)Hi " 6!d 2^ ( _ a) .i, " 6!d aFl( " 6 ' " C ' " a ' 1} ' 

i=0 v ; v ; i=0 v ' 

where we use that the sum is an 2 Pi evaluated at 1. Its value is given by the 
Chu-Vandermonde formula, see pQ, hence 

1V - (a-»)! _ a!(c- o) 6 
f^ 1 j (6-*)!(c-t)!i! (-o) 6 6!c! ' 

The two sums above are of this form and we get 

[n — k — r)\(k + 1 — rj! 



where 



Q=(2fc + 2) (2fc ~ 2n + 1 ^- fc --(fc + l-r) ( 2fc - 2 ^- 



(fc + r - 2n) n _ fc _ r (fc + r - 2n) n _ fc _ r 

= (2k -2n + l)n-k-r-i r, 2k + 2){k _ r _ n) _ {k + l _ r w 2k _ 2n)] 
(k + r - 2n.) n _ /c _ r 

= 2r(ra + 1) , 

where we used (a) n = (— l) n (l — a — n) n twice. This gives 

n + 1 \(2n-2Jfe-l)! 



(-l)V r) (l)=2r^ 



Jc+l — rJ (n — k — r)\ 
and finally 

(n-fc)A(fc+l) . . , , 

/- , r 2r/ n + 1 \ 2n-2fc-l l 

PW= > I-" 1 , , 1 )—( Z \T- 

r\\k + \ — rl (n — k — r)\ 

Note that the term corresponding to r = is zero. If we insert this expression 
for p(a) in (|19|1. we get the formula of Theorem 12.11 □ 

Remark 3.3 The evaluation above of (— l) r p^(l) can be done using generat- 
ing functions like in The authors want to thank Mogens Esrom Larsen for 
the idea to use the Chu-Vandermonde identity twice. 
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3.3 Proof of Theorem EO 

The starting point is the following formula of Macdonald, see ^2] 

s z 2 + X 2 zX ds 

— \K V { ) — , 

2s s s 



1 f°° s 
K u {z)K v {X) = -J exp[-~ 



(20) 



which we will use for v = n + i, z = au, X = (1 — a)u. Multiplying ()2()fi by 



' (a(l-o)« 2 )" 



and using © we find 
1 



2T(i/) Jo 



exp 



-2~ U 



I> 

k v (au)k u ((l — a)u 
,a 2 + (l-a) 2 



2s 



s h v 



a(l — a)u z 



ds. 



We now insert that with v = n + \ we have = e '"'(?n(|^|) and hence 

after some simplification 

e-\ u \q n {a\u\)q n ({l-a)\u\) = 



ds. 



2 "+!r(n + i) Jo 



poo 


s 


u 2 - 




/ exp 
Jo 


~2 


2s _ 


--^( a(1 - a)u2 ) 



We next insert the expression Q for q n under the integral sign. This gives 

e~^q n (a\u\)q n ((l - a)\u\) = 
1 



J>fV(l-a)) 



k=0 " _i V" i 2' 

Using the following formula from 7, 3.471(9)] 



poo 


s 


u 2 " 


/ exp 
/o 


~2 


2s 



jn—k- 



2 ds. 



X 



x u ~ l exp(-^ - 7 a;) dx = 2 ( £ ) (2^ 
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(21) 



and again © the above is equal to 

.2- k +h(n-k + l) -\u 



J>f(a(l-a))* 



fc=0 



2™+2r(n + |) 



< '-l« 2fc g n _ fc (|u|). 



Finally, using Pochhammer symbols and skipping absolute values since we are 
now dealing with a polynomial identity, we get 



: ( 9 k 



q n (au)q n ((l - a)u) = Ja[ B, (a(l - a)) fc ^V fc g n _ fe ( U ). (22) 



fc=0 
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Using the expression for u 2k q n -k{u) from Lemma 13. II and the expression for 
Ofc in (|23) we then get 

q n (au)q n ((l - a)u) = E S ui ~ °))* E^W") 



i=o fc=i UM2^ ft ' 



hence 



with 



g n (au)<? n ((l - a)u) = ^ j3\ n) {a)q n+i { 

i=0 



k 

n—i 



1, , 1, 



Collecting 
we get 



x ( n -i-Z+-), +aj (-n- -)/ 



(-) n _i_j(n - i - J + -)/ +ai = (-)«+< 



# » = (a(l - a)) ^(a(l - a)) 2 TTTT^ " 9)' 



wl _ o)) .^%^£ (4o(1 _ o) yP"--0K-»-J), 



n! 

(2n)!(i) n t!^""^ ~ y/ (n-i-Z)!Z! 



n! \ 2 (2n + 2i)!^, 2 (2n - i - 1)1 (-n - £), 



v v V(2n)! / + _ 



Expanding (1 — (2a — l) 2 )' using the binomial formula and interchanging the 
sums we get 

1=0 v ; 

^ • (2a-lg> ^ (2n-»-Q!(-n-|), 

^ j j! t; (n-i-iy.(i-j)i ■ 
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We claim that 

' (2n-i-l)l(-n-h)i 



s> = E 



- Iff ^(-«K 2 " 2 r)(T')- <»> 

and we have obtained the final formula 

#»(<,) = (mi - (^y 2 - ^ 2 "- 2i >'< 2 " +2i > ! 



(2n)\J (n-i)\(n + i)\ 



«g(% +1 )("7>-^ 

We will see that (j2Hjl is a Chu-Vandermonde formula. In fact, shifting the 
summation index putting / = j + m we get 



1 \V(2n-i-j-m)!(-n + j-±) 



2 ;m 



2 m!(n — « — j — m)\ 

and calling the general term in this sum c m we get 

c m +i _ {m + i+ j -n)(m + j -n-\) 



c m (m + l)(m + i + j — 2n) 

which shows that the sum is a 2-Fi- We have 

S = (-n- \)jC Q x 2Fi(-(n-i- j),j -n - -;i + j - 2n;l), 

and using the Chu-Vandermonde identity, cf. pQ: 

(c - a) n 



iFi(-n,a;c; 1) 



(c)n 



we get 



1 (2n-i- j) \ (-n + z + \) n -i-j 

2 J (n-i-j)l (i + j - 2n) n _ i _ i 

_ Is (2n-i-j)!(j + ^)n-i-j 
' " 2 Jj (n-i- j)! (n + l) n . 



where we used {a) n = (— l) n (l — a — n) n twice. We can now simplify to get 
q _ ( _ 1 \ n] - (gV-i 

and applying the formula 
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twice we get 

(2n-2i)! 2i _ 2 _ 1 n! j!2 2 ^ 

(n-i)! 1 2 ;j (n-i- j)! (2j)!' 

Now we can write 

( n 2 jj2 " ( ij n!(2n-2j + l)!' 

and hence 

(n-i)\ V 2 J 

□ 



3.4 Proof of Theorem EOl 

For n, m > and o £ i, we can write 

m+n 



g n (au) q m ((1 - a) «) = ^ /3^' m) (a) % (n) (24) 

fe=0 



for some uniquely determined coefficients since the left-hand side is a polynomial 
in u of degree <n + m. Clearly /?[ n '"^ (a) is a polynomial in a satisfying 

^ m \a)=(3l m ' n) (l-a). (25) 

We shall prove that (3^' m) (a) > for < a < 1 and that /?[ n,m) (a) = if 
k < n A m, which will be a consequence of the following recursion formula. 

Lemma 3.4 For n,m> 1, we have 



1 ,q("' m ) _ fl2 a(n-l,m) , n. (1 «) a (n,m-l) 

2A , + 1 ^ +1 W- 2n _!^ W+ 2m-l^ 



where k = 0, 1, . . . ,m + n — 1. 
Furthermore Pq 1 '™^ (a) = 0. 

Proof: Differentiating (|24|) with respect to u gives 

m+n 



a<?„ (an) <? m ((1 -a)u) + (l- a) q n (au) q' m ((1 -a)u)= ^ /3^™' m) (a) q£ (u) 



,(n,m) ^ ^ 

fc=l 

and using the formula (|14|) we find 

a (a„ (au) - (a«)j <? m ((1 - a) u) 

+ (1 - a) g n (an) fg m ((1 - a) u) - ((1 - a) «)J 
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and using (|24[) once more we get 



-^TiQn-i {au) q m ((1 - a)u) - ( 2r ^\ u q n (au) q m ~\ ((1 - a) u) 

n{n,m) I \ ,-^ n -|_ m _i a (n,m) , \ /„, in— 1 i \ 

-0o ( a )- u 22k=o Pk+i '(a)(2fc + l) Qk(u). 

For u = this gives ft^ 1 (a) = and dividing by — u and equating the coeffi- 
cients of qy. (u) , we get the desired formula. □ 

Now the proof of Theorem l2,4l is easy by induction in k and by the symmetry 
formula (|25j). we can assume n > m. Let < a < 1. We prove that/4"' m) (a) >0 
for k < n + m and that it is zero for k < m (under the assumption n > m). 
This is true for k = by Lemma 13.41 when m > 1, and for m = it follows 
by Theorem 12.11 Assume now that the results hold for k = k$ and assume 
ko + 1 < n + m. The nonnegativity for A; = fco + 1 now follows by Lemma l3.4[ 
and likewise if ko + 1 < m < n the coefficient is since ko < (n — 1) A (m— 1). 
□ 

3.5 Proof of Theorem E31 

By Theorem 12.41 the result holds for k = 2. Assuming it holds for k — 1 > 2 we 
have 

L' 

9m ' ' ' ?n fc _i (afc_iu) = ^7i9i((l -a>k)u), ueR (27) 

with /' = min(ni, . . . , nk~i),L' = m + • • • + Wjt-i and 7,- > because we can 
write 

ajU=- — - — (l-ajfc)u, i = 1, . . . ,k - 1. 
1 — dk 

If we multiply (|18JI with q nk {aku) we get 

L' V n k +j 

^2ljQn k (aku)qj((l - a k )u) = ^7j X] &i (. a k)<li( u ), 
3=1' 3=1' i=n k /\j 

and the assertion follows. □ 

4 Inverse Gamma distribution 

Grosswald proved 8 that the Student-t distribution is infinitely divisible. This 
is a consequence of the infinite divisibility of the inverse Gamma distribution 
because of subordination. It was proved later that the inverse Gamma dis- 
tribution is a generalized Gamma convolution in the sense of Thorin, which 
is stronger than self-decomposability and in particular stronger than infinite 
divisibility. 

The following density on the half-line is an inverse Gamma density with 
scale parameter | and shape parameter u > 0: 

C,exp(-i)t— \ C u = ¥ ^— (28) 
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Let the corresponding probability measure be denoted % and let 

1 x 2 
V47rt 44 

denote the Gaussian semigroup of normal densities (in the normalization of [5] ) . 
Then 

f v (x) = / g t (x)d%{t) (29) 



is the Studcnt-t density © with 2v degrees of freedom. The corresponding 
probability measure is denoted a v . This formula says that a u is subordinated 
to the Gaussian semigroup by an inverse Gamma distribution, and it implies the 
infinite divisibility of Student-t from the infinite divisibility of inverse Gamma. 
Since the Laplace transformation is one-to-one, it is clear that if two probabil- 
ities 71,72 on ]0,oo[ lead to the same subordinated density 

g t (x)dji(t)= / gt(x)dj 2 (t), xeM., 
Jo Jo 

then 71 =72. 

If we denote T a (x) = ax, the distribution r a (a 1 ) * Ti_ a (<r 1 ) is given in 



2 



m+7 



(fTU|) . However note that T a (g t (x)dx) = g ta i{x)dx so 

r a ((Tv)= g ta 2(x)d%(t) dx, (30) 
J 

hence 



T a (cr Vl ) * Ti_ a (o^) = / / (g ta 2dx) * (# s (i_ a )2Gfe) d%i {t)dju 2 ( s ) 
Jo Jo 



00 roo 



JO 

00 



{g t a^+s(l-a)2dx) d% x (t)d% 2 (s) 



g u {x) dr a 2^ 

2 (7^)(«) dx. 

Therefore, using (|3U|) we see that for i>i = n + ^,t / 2 = m+ ^ with n, m = 0, 1, 
the formula H1U|) rewritten as 

n+m 

r a (cr 1) *ri_ a (fj 1) = S~] W (a)o - , , 1 

n+2 m +2 ^""^ + 2 

k=n/\m 



is equivalent to 

n+m 

a (n,m) / n~ 

0* ( a)7 <--2 



^ 2 (7 n+ i)*T(i-a) 2 (7 m+ i)= E 4 n ' m) («)7 fc+ i- (31) 



k=nhm 



This shows that Theorem 12.41 is equivalent to the following result about inverse 
Gamma distributions: 



14 



The distribution of a 2 Z n + (l — a) 2 Z m , where Z n , Z m are independent inverse 
Gamma random variables with distribution \2ty) for u = n+^,m+^ respectively, 
has a density which is a convex combination of inverse Gamma densities. 

This result can be extended to the multivariate Student-t distributions as 
follows. A rotation invariant iV— variate Student-t probability density is given 
for x = (xi, . . . , xn) 6 M. n by 



f N , u (x)=A N , u {l+\x\ 2 ) 2 , A N>U = - 7 -l2Z_ 



rM (r(D) 

where 

N I 
(x,y>=^x i y i , |x| = ((x,x))2 , x,yel ff . 

i=i 

It is easy to verify that /^^(x) is subordinated to the N- variate Gaussian 
semigroup 

TV / | X | 2 \ 

g N}t (x) = (47rt)~2 exp f ~^-\ , t>0,x€R N 
by the inverse Gamma density ()28|) . i.e. 

poo 

f N ,u(x)= / 9N,M)d7 v (t). (32) 
Jo 

Therefore the characteristic function is given by 



/ e J < x ^/ JV ,,(x)dx = ^(|y|) (33) 

generalizing Q. In fact 

e l{x ' y) f N ,„ (x) dx = f°° ( [ e^g N , t (x) dx) d%{t) 
Jo \Jm N J 



oo 

t\y\ 2 



o 



and the result follows by (|2"T|) . 

As a conclusion, the Theorems 2.1, 2.2 and 2.4 apply in the multivariate 
case. For example, an equivalent form of (|1U|) writes as follows: with < a < 1, 

^ ^ n+m 

^v//v,n+i (a~ lx )* ! yv -ftv.m+i (C 1 - x ) = X! ^"' m) /jv,fc+| ( x ) 

\ ' k=nAm 
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